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ABSTRACT 

The combined effect of a transverse magnetic field and radiative heat transfer on unsteady flow of a 
conducting optically thin visco-elastic fluid through a rotating parallel channel filled with saturated porous 
medium and non-uniform walls temperature has been discussed. It is assumed that the fluid has small electrical 
conductivity and the electromagnetic force produced is very small. The analytical solutions are obtained for the 
problem making use of perturbation technique. The effects of the radiation and the magnetic field parameters on 
velocity profile and shear stress for different values of the visco-elastic parameter with the combination of the 
other flow parameters are illustrated graphically, and physical aspects of the problem are discussed 
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I. INTRODUCTION 

In recent years, there has been a considerable interest in rotating hydro magnetic fluid flows due to 
possible applications to geophysical and astrophysical problems. The magnitude analysis shows that in the 
basic field equations, the Coriolis force is very significant as compared to the inertial force. Furthermore, it 
reveals that the Coriolis and magneto hydro dynamic forces are of comparable magnitude. It is generally 
admitted that a number of astronomical bodies (e.g. the Sun. Earth. Jupiter, Magnetic Stars, pulsars) possess 
fluid interiors and (at least surface) magnetic fields. Changes in the rotation rate of such objects suggest the 
possible importance of hydro magnetic spin-up. The flow of an electrically conducting fluid has important 
applications in many branches of engineering science such as magneto hydro dynamics (MHD) generators, 
plasma studies, nuclear reactor, geothermal energy extraction, electromagnetic propulsion, and the boundary 
layer control in the field of aerodynamics. In the light of these applications, MHD flow in a channel has been 
studied by many authors; some of them are Nigam and Singh [1], Soundalgekar and Bhat [2], Vajravelu [3], and 
Attia and Kotb [4]. A survey of MHD studies in the technological fields can be found in Moreau [5]. The flow 
of fluids through porous media is an important topic because of the recovery of crude oil from the pores of the 
reservoir rocks; in this case, Darcy's law represents the gross effect. Raptis et al. [6] have analysed the hydro 
magnetic free convection flow through a porous medium between two parallel plates. Aldoss et al. [7] have 
studied mixed convection flow from a vertical plate embedded in a porous medium in the presence of a 
magnetic field. Makinde and Mhone [8] have considered heat transfer to MHD oscillatory flow in a channel 
filled with porous medium. Recently the combined effect of a transverse magnetic field and radiative heat 
transfer on unsteady flow of a conducting optically thin visco-elastic fluid through a channel filled with 
saturated porous medium and non-uniform walls temperature has been discussed by Rita Choudary and Utpal 
Jyothi Das [12]. In this study, an attempt has been made to extend the problem studied by Makinde and Mhone 
[8] to the case of visco-elastic fluid characterised by second-order fluid. The combined effect of a transverse 
magnetic field and radiative heat transfer on unsteady flow of a conducting optically thin visco-elastic fluid 
through a rotating parallel channel filled with saturated porous medium and non-uniform walls temperature has 
been discussed. It is assumed that the fluid has small electrical conductivity and the electromagnetic force 
produced is very small. The analytical solutions are obtained for the problem making use of perturbation 
technique. 
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II. MATHEMATICAL FORMULATION AND SOLUTION OF THE PROBLEM 

Consider the flow of a conducting optically thin fluid in a channel filled with saturated porous medium 
under the influence of an externally applied homogeneous magnetic field and radiative heat transfer. It is 
assumed that the fluid has small electrical conductivity and the electromagnetic force produced is very small. 
The x-axis is taken along the centre of the channel, and the z-axis is taken normal to it. In the initial undisturbed 

state both the plates and the fluid rotate with the same angular velocity O. . At t > 0 , the fluid is driven by a 
constant pressure gradient parallel to the channel walls. Then, assuming a Boussinesq incompressible fluid 
model, the equations governing the motion are given by 

du I dp 8 2 u 8' u a B 2 u u 

2Q v = + V| — + v 2 — r -- •- -v, — + g/3(T -T 0 ) (1) 

8t p 8x 8z 8z 8t p K 

2 3 2 

8v I dp d v 8 v a B„v v 

+2nu = - +v, + v, - 0 - v , — (2) 

12 2 2 1 

8t p 8y dz 8z 8t p K 

dT k d 2 T 1 dq 

— = r - (3) 

8t pC dz pC dz 
Corresponding boundary conditions 

u = 0,v = 0,7 = T w on z = 1 (4) 

u = 0,v = 0,T = T o on z = 0 (5) 

Where u is the axial velocity, t is the time, T is the fluid temperature, P is the pressure, g is the gravitational 
force, q 1 is the radiative heat flux, /? is the coefficient of volume expansion due to temperature, C is the 
specific heat at constant pressure, k is the thermal conductivity, 

K is the porous medium permeability co-efficient, B o (= p H o ) is electromagnetic induction, fj is the 
magnetic permeability, H is the intensity of the magnetic field, <r is the conductivity of the fluid, p is fluid 

density, and v . = p . I p , (i = 1, 2) . It is assumed that both walls of temperature T g , T are high enough to 

induce radiative heat transfer. Following Cogley et.al [11], it is assumed that the fluid is optically thin with a 
relatively low density and the radiative heat flux is given by 

dq , 

'-= 4a;(T o -T), (6) 

dz 

Where a , is the mean radiation absorption co-efficient. 

Combining equations (1) and (2) and let q = u + iv and <f = x + iy , we obtain 

dq 1 dp d 2 q d^q a B^q q ,„ m 

-^-+2iClq = l -+v l —^+v 1 ^--^^--v l ^-+ gf}(T-T 0 ) (7) The 

dt p d% dz dz dt p K 

following non-dimensional quantities are introduced: 

x y , £, , z u , v 

x = — , y = — , £ = — , z = — , u = — , v = — , 
a a a a U U 

1 tU , ap T - T 0 

q = — , t = — , p = , & = 

U a pv p T w - T 0 

Making use of non-dimensional variables, the dimensionless governing equations together with appropriate 
boundary conditions (dropping asterisks) are 

dq dp 8 2 q S 3 q , _, , 

Re = + -+ a (M " + HE + S )q + Gr T (8) 

dt 8% dz 2 8z 2 8t 

80 d 2 0 

Pe = + R 0 (9) 

dt dz 1 
With 
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q = 0, 0 = 1 

q = 0, 0=0 
where 

Ua 



on 
on 



z = 1 

z = 0 



(10) 
(11) 



Re 



is the Reynolds number, 



M 



a B a a~ v u 

— is the Hartmann number, E = — is the Eckmann number 



v.U 



Cla' 



k 1 
D = — is the Darcy number (or) S = — is the porous medium shape factor parameter, 

a 2 D 

v Re 

a = — is the visco-elastic parameter, 



Gr 



g/3(T w -T 0 )a 2 . 



is the Grashoff number, 



Ua pC 

Pe = is Peclet parameter and 



R 



k 

4a , a ' 



is the Radiation parameter. 



solving the equations (2.8) and (2.9) for purely oscillatory flow, Let 

dP 



Ae 



dx 



q(z,t) = q 0 {z)e'°" , 0(z,t) = 0 o (z)e"" 



(12) 



where, A. is constant and w is the frequency of oscillation. 

Substituting the above expressions (12) into the equations (8) and (9), and making use of the corresponding 
boundary conditions (10) and (1 1), we obtain 



(1 + iaw ) m t q 

dz 1 



-Z - Gr 



— + m;0 o = 0 



d 0. 

dz'' 

subjected to the boundary conditions 
o„ = 0,0=1 on z = 1 



(13) 
(14) 



(15) 

q o = 0, 0 O = 0 on z = 0 (16) 

where m 1 = yj M 2 + 2iE _1 + S 2 + ico Re and m 2 = yj R 2 - ico Pe 

Equations (13) and (14) are solved; we obtained the solution for the fluid velocity and temperature as follows 



q{z,t) = \ a, e"" 1 "'"' 
I 

sin( m.z) 

0( z ,t) = 5— e ' 

sin( m ) 



A 



A 



Gr sin( m 2 z ) 
(m 2 b + m ' ) sin( m , ) 



(17) 
(18) 



where a 



< A > 



Gr 

(m 2 b + m 2 ) 



- , b = 1 + i am 
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The non-dimensional shear stress a at the wall z = 0 is given by 

a ( du d u 

I — + 



a 



{up la) {dz SzdtJ^ o 



T a in 
I — 

L b 



Gr 



1 

| ( 1 + ( ooa ) e ' 



(m t b + m l ) sin( m 2 ) J 
The rate of heat transfer across the channel wall z = 1 is given as 

f & ^ 1 m , C0S ( m , ) tml 



Nil 



sin( m ) 



(19) 



(20) 



III. RESULTS AND DISCUSSION: 

The flow governed by the non dimensional parameters, Re is the Reynolds number, M 2 is the 
Hartmann number, E Eckmann number, D is the Darcy number (or) S is the porous medium shape factor 
parameter, a is the visco-elastic parameter, R is the Radiation parameter with fixed values of Gr the 
Grashoff number, Pe Peclet parameter. We have considered the real and imaginary parts of the results u and v 
throughout for numerical validation. The velocity profiles for the components against z is plotted in Figures 1— 
12 while figure 13-18 to observe temperature profiles on the visco-elastic effects and other parameters for 
various sets of values of Hartmann number H, porous parameter S and radiation parameter R, with fixed values 
of other flow parameters, namely, Pe = 2, t = 0.1, Gr = 2, 1 = 1, and co = 1. 

It is evident from Figures (1-12) that the velocity profiles is parabolic in nature, and the magnitude of 
velocity u and v increase with the increasing values of the Reynolds number Re, Porous parameter S, the visco- 
elastic parameter I a I , Radiation parameter R and (Figure 1, 2, 5-12). It is also noted from the figures (3-4) 
that the magnitude of the velocity component u experiences retardation and the behaviours of the velocity 
component v remains the same with the increasing values of the Hartmann number. We observe that lower the 
permeability of the porous medium lesser the fluid speed in the entire fluid region. The resultant velocity q 
enhances with increasing the parameters Re, D,a , R and experiences retardation with increasing the intensity of 
the magnetic field. It is evident that the temperature profiles exhibit the nature of the flow on governing 
parameters. The magnitude of the temperature increases with increasing Re, D,a , R, E and experiences 
retardation with increasing the magnetic field parameter (Hartmann number M) (Figures 13-18). 

The shear stress on the wall and the rate of heat transfer evaluated analytically and computationally 
discussed with reference to governing parameters (Tables 1-2). It is evident that the shear stress enhances with 
increasing Re, M and a and reduces with increasing S. We also noted that it is increases firstly and then 
decreases with the increasing values of radiation parameter R. Table. 2 depict that the Nusselt number (Nu) 
increases with the increasing values of the visco-elastic parameter I a I and the magnitude of the heat transfer 
increases with increasing Re, S, a , R, E retardation with increasing the magnetic field parameter. It has also 
been observed that the temperature field is not significantly affected by the visco-elastic parameter. 
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Fig. 1 The velocity profile for u on Reynolds number 

Re 

with 



Fig. 2: The velocity profile for v on Reynolds number 
Re with 
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-0.1, £ = 0.01,/? = 1.5 
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Fig. 3: The velocity profile for u on Hartmann number 
M with 

Re =50, S =1, a = -0.1, E = 0.01,2? = 1.5 



Fig. 4: The velocity profile for v on Hartmann number 
M with 



Re = 50 , S = 1, a 



-0.1, E = 0.01, R = 1.5 




1 

0.5 
0 




0.5 
z 



-♦— S=l 
-■-5=2 

-*-S=3 

S=4 



Fig. 5: The velocity profile for u on Porous parameter 
S with 

Re = 50, M = 2, a = -0.1, E = 0.01,/? = 1.5 



Fig. 6: The velocity profile for v on Porous parameter 
S with 

Re = 50 , M =2, a = -0.1, E = 0.01,7? = 1.5 
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Fig. 7: The velocity profile for u on visco-elastic 
parameter a with 

Re =50, M = 2, S = I, E = 0.01,7? = 1.5 



Fig. 8: The velocity profile for v on visco-elastic 
parameter a with 

Re =50, M =2, S = l,E = 0.01,7? = 1.5 
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Fig. 09: The velocity profile for u on Radiation 
parameter R with 

Re =50, M =2, S = 1, E = 0.01, a = -0.1 



Fig. 10: The velocity profile for v on Radiation 
parameter R with 

Re = 50 , M =2, S = 1, a = -0.1 , E = 0.01 
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Fig. 1 1 : The velocity profile for u on Eckmann number 
E with 

Re = 50, M = 2, S = 1, E = 0.01, R = 1.5, a = 



Fig. 12: The velocity profile for v on Eckmann 
number E with 

Re =50, M =2, S = 1, E = 0.01, R = 1.5, a 
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Fig. 13: The temperature profile for 9 on Reynolds 
number Re with 



M = 2, S = 1, a 



-0.1, £ = 0.01, R = 1.5 



Fig. 14: The temperature profile for 9 on Hartmann 
number M with 

Re = 50, 5=1, a = -0.1, £ = 0.01, R = 1 .5 
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Fig. 15: The temperature profile for 9 on Porous 
parameter S with 
Re =50, M =2, a = -0.1, E = 0.01 , R = 1 .5 



Fig. 16: The temperature profile for 9 on visco- 
elastic parameter a 
with 

Re =50, M = 2, 5=1, £ = 0.01,tf = 1.5 
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Fig. 17: The temperature profile for 


9 on Radiation 


Fig. 18: The temperature profile for 8 on 


Eckmann 


parameter R 


with 


number E 


with 


Re =50, M = 2, S = 1, a = - 


0.1, £ = 0.01 


Re =50, M = 2, S = 1, a = -0.1 , E 


= 0.01, R = 



Re 


I 


II 


III 


IV 


V 


VI 


VII 


VIII 


IX 


X 


XI 


50 


0.0245 


1.8554 


3.8809 


0.0126 


0.0052 


0.0546 


0.0785 


0.0854 


0.0569 


0.1055 


0.2551 


100 


0.1446 


2.8847 


5.8824 


0.1222 


0.0145 


0.1545 


0.1884 


0.1885 


0.1254 


0.2526 


0.5887 


150 


0.2114 


5.9968 


8.1055 


0.1884 


0.1225 


0.2405 


0.2887 


0.2995 


0.2254 


0.5855 


0.7886 


200 


0.2889 


9.5548 


12.658 


0.2254 


0.1889 


0.2998 


0.3225 


0.3825 


0.2996 


0.6335 


0.8888 





I 


II 


III 


IV 


V 


VI 


VII 


VIII 


IX 


X 


XI 


M 


2 


5 


8 


2 


2 


2 


2 


2 


2 


2 


2 


S 


1 


1 


1 


2 


3 


1 


1 


1 


1 


1 


1 


a 


-0.1 


-0.1 


-0.1 


-0.1 


-0.1 


-0.2 


-0.3 


-0.1 


-0.1 


-0.1 


-0.1 


R 


1.5 


1.5 


1.5 


1.5 


1.5 


1.5 


1.5 


2.5 


3.5 


1.5 


1.5 


E 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.02 


0.03 



Table. 1 The shear stresses ( r ) at the wall z = 0 



Re 


I 


II 


III 


IV 


V 


VI 


VII 


VIII 


IX 


X 


XI 


50 


1.2552 


0.9885 


0.6585 


2.3314 


3.9985 


2.8845 


3.1145 


2.0147 


3.0054 


1.8225 


1.9969 


100 


2.0458 


1.5665 


1.1125 


3.2214 


4.2256 


3.2251 


4.6556 


3.1022 


4.0552 


2.2258 


2.5566 


150 


2.9985 


2.4458 


1.8447 


3.9985 


6.5582 


4.8859 


5.9985 


3.5229 


4.8002 


3.2556 


2.2985 


200 


3.6687 


3.1205 


2.7195 


4.8556 


7.5222 


5.6648 


6.7808 


4.2245 


5.1566 


3.8853 


4.1205 
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II 


III 


IV 


V 


VI 


VII 


VIII 


IX 


X 


XI 


M 


2 


5 


8 


2 


2 


2 


2 


2 


2 


2 


2 


S 


1 


1 


1 


2 


3 


1 


1 


1 


1 


1 


1 


a 


-0.1 


-0.1 


-0.1 


-0.1 


-0.1 


-0.2 


-0.3 


-0.1 


-0.1 


-0.1 


-0.1 


R 


1.5 


1.5 


1.5 


1.5 


1.5 


1.5 


1.5 


2.5 


3.5 


1.5 


1.5 


E 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.02 


0.03 



Table. 2 The Rate of heat transfer (Nusselt number) at the wall z = 1 



IV. CONCLUSIONS: 

1. The magnitude of velocity u and v increase with the increasing values of the Reynolds number Re, 
Porous parameter S, the visco-elastic parameter a , Radiation parameter R and Eckmann number E. 

2. The magnitude of the velocity component u experiences retardation and the behaviours of the velocity 
component v remains the same with the increasing values of the Hartmann number. 

3. Lower the permeability of the porous medium lesser the fluid speed in the entire fluid region. 

4. The resultant velocity q enhances with increasing the parameters Re, D, a , R, E and experiences 
retardation with increasing the intensity of the magnetic field. 

5. The magnitude of the temperature increases with increasing Re, D, a , R, E and experiences retardation 
with increasing the magnetic field parameter (Hartmann number M) 

6. The shear stress enhances with increasing Re, M and a and reduces with increasing S. We also noted 
that it is increases firstly and then decreases with the increasing values of radiation parameter R. 

7. The Nusselt number increases with the increasing values of the visco-elastic parameter I a I and the 
magnitude of the heat transfer increases with increasing Re, S, R, E retardation with increasing the 
magnetic field parameter M. 
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